a^ 
a^ 

CD 

m 

(N 

(N 
> 

o 

o 

a^ 



X 



Non-perturbative improvement of bilinears in unquenched QCD * 
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We describe how the improvement of quark bilinears generalizes from quenched to unquenched QCD, and 
discuss which of the additional improvement constants can be determined using Ward Identities. 



A major motivation for undertaking the im- 
provement program is to facilitate unquenched 
calculations by allowing simulations at larger lat- 
tice spacings. Thus it is important to study the 
application of the improvement program to un- 
quenched theories. This has been done for the 
action itself, and for certain operators in the chi- 
ral limit ||l|,||. Here we take a further step by 
considering the theory of on-shell improvement 
of quark bilinears in unquenched QCD with non- 
zero quark masses. We enumerate the additional 
improvement coefficients that are required, and 
discuss which of them can be determined non- 
perturbatively using Ward Identities (WI). 

The analysis depends on the number of dynam- 
ical flavors, and we consider here the physically 
relevant theories with Nf > 3 non-degenerate fla- 
vors. The two flavor theory is more complicated 
and will be discussed elsewhere We use an 
abbreviated notation for flavor traces 



Aijj) 



(AA„) = Tr(AA^), 



with A an SU{Nf) generator. Note that we con- 
sider both flavor singlets and non-singlets, with 
the latter being both off-diagonal and diagonal. 
All three are needed for phenomenology. For ex- 

('23) 

ample, the flavor off-diagonal operator ' — 
d'jfj.jss is needed to determine fx, and the nu- 
cleon matrix elements of the flavor diagonal op- 
erators Ajf and T^^ (which, with j an unsummed 
flavor index, are linear combinations of flavor 
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singlet and non-singlet) give information on the 
structure functions. 

We begin by reviewing previous work on non- 
perturbative 0(a) improvement of unquenched 
QCD. The ALPHA collaboration g] has shown 
how on-shell improvement of the action can 
be accomplished by adding the Sheikholeslami- 
Wohlert or "clover" term, with appropriately cho- 
sen coefficient csw- On-shell matrix elements of 
flavor non-singlet axial, vector and tensor bilin- 
ears are then improved, in the chiral limit, by 
adding dimension four operators: 

{XA^y - {XA^)+CAd^{XP), (1) 
{XV^y = {XV^)+cvd,{XT^,), (2) 
(AV)^ = (AV>+CT[a^(AK)-a,(AF^)].(3) 

Here we use the standard notation for local lattice 
bilinears (see |jl|,Q), and, for brevity, set the lat- 
tice spacing a to unity. Flavor non-singlet scalar 
and pseudoscalar bilinears are automatically on- 
shell 0{a) improved in the chiral limit. 

The improvement coefficients csw and cp (F = 
V, A, T) depend on the number of dynamical 
quarks, Nf, and on the bare coupling g^. They 
can all be determined non-perturbatively using 
WI. In particular, enforcing the partial conser- 
vation of the improved flavor off-diagonal axial 
current determines csw and ca, as well as the 
critical quark mass nic at which chiral symme- 
try is restored |l| . This method works for two or 
more flavors, and has been implemented numer- 
ically for Nf = 2 With the improved axial 
current in hand, one can then enforce the axial 
transformation properties of the vector and ten- 



sor bilinears to determine cy |5|J^ and ct ||§1- 

Away from the chiral limit there are many ad- 
ditional improvement coefficients. In particular, 
the effective coupling constant becomes ||] 



l^=9l{l + b,{M)/Nf), 



(4) 



where (M) is the trace of the mass matrix. 
Knowledge of the improvement coefficient bg al- 
lows one to adjusts go, as (M) is varied, such that 
5off is constant. This corresponds to working at 
fixed lattice spacing. In this way one does not 
introduce spurious 0{aM) dependence in physi- 
cal quantities. Methods for determining bg non- 
perturbatively are given in Refs. [0J^. 

This completes our review of previous work. 
We now describe the new improvement coeffi- 
cients that arise when one improves bilinears in 
the unquenched theory. First we note that, in 
the chiral limit, improvement of flavor singlet bi- 
linears is more complicated than that of non- 
singlets, due to the contributions from discon- 
nected contractions. For axial, vector and tensor 
bilinears the form remains the same, e.g. 



(5) 



with new constants ca, cy and ct- For the other 
bilinears there are additional gluonic operators^ 



(6) 
(7) 



In principle, one can determine all five new im- 
provement coefficients by enforcing the invariance 
of these operators under an off-diagonal vector 
transformation, i.e. Sy'^^O^ = 0. This equation 
can be made non-trivial by taking appropriate fla- 
vor non-singlet matrix elements. 

Moving away from the chiral limit, there are 
several new improvement terms. For flavor non- 
singlets, the general form is exemplified by 



Za 



{1 + bA{M))(XA, 



+bA^{{X,M}A^) + fA{XM){A^)] , 



(8) 



where the "hat" indicates an improved and renor- 
malized operator. CP invariance implies that 

'Normalized such that 5giue = (l/2gQ)tr(_F^yi<)j,y). 



only the anticommutator of A and M appears. 
It is instructive to consider examples of this gen- 
eral formula. For flavor off-diagonal operators, 
the /-term drops out, leaving 

4^") = [l + bA{M) + bA^im,+m2)] 

The bA term is that present in the quenched ap- 
proximation, while bA multiplies the additional 
mass dependence arising from quark loops. All 
terms contribute for diagonal operators, e.g. 



A 



(11) 



-A 



(22) 



(4")-4^^))^(l + 6^(Af)) 



+bA{nHAl''^~m2A<'^^^) + fA{mi-m2){A^) 



The /-term arises from disconnected contractions 
of the operator, and is present only for non- 
degenerate quarks. 

Finally, we consider the mass-dependent im- 
provement coefficients needed for flavor singlet 
operators. Here there are only two independent 
traces, and thus two terms per bilinear, e.g. 

(A^) = Zata [(l + dA(M))(A^)VdA(AfA^)] (9) 

Note that the overall normalization constant dif- 
fers from Za by a factor ta- For the axial current, 
this factor is scale-dependent since the singlet cur- 
rent has a non-zero anomalous dimension. For 
the other bilinears, however, rr is a finite, scale- 
independent quantity, dependent only on the ef- 
fective coupling constant and on N f. This is be- 
cause the anomalous dimensions of the singlet and 
non-singlet bilinears are the same. 

In summary, the number of improvement terms 
involving mass dependence increases substan- 
tially when one unquenches and considers both 
singlet and non-singlet bilinears. Improvement 
of flavor off-diagonal operators in the quenched 
theory requires the determination of 5 such co- 
efficients, while complete improvement of bilin- 
ears in the unquenched theory requires 24. In the 
quenched theory, all 5 coefficients can be deter- 
mined using WI with non-degenerate quarks 
We have studied the generalization of this method 
to unquenched QCD, and find that all but three 
of the new coefficients can be determined using 



WI with the following steps. ^ 

1. Enforcing the conservation of the diagonal vec- 
tor charges while independently varying the three 
quark masses determines all five coefficients asso- 
ciated with the vector bilinear (as well as ry). 

2. Enforcing the PCAC relation away from the 
chiral limit (combined with considerations out- 
lined below) determines the combinations bs, 
bp-bA, bp-bA-bs, fs, ds-^bs, and ds-bs, as 
well as rs ■ This generalizes the method of Ref. . 

3. Enforcing the off-diagonal axial transformation 
properties of operators, for non-degenerate quark 
masses. This is the generalization of the method 
of Ref. iQ, and determines by + bA,bv + bA, bp + 
bs, '2^/rsibp - bs) + bp - bs, and bx- 

4. Enforcing off-diagonal vector transformation 
properties of operators, e.g. 

412)0(21) ^ 0(11) „ 0(22)^ Si^^)(^0)^0. (10) 

These identities have not been considered pre- 
viously, since they involve flavor diagonal bilin- 
ears. They determine the /r (T = V, A, T, S, P) 
together with cLa, dp and cLt, and rp. 

5. Enforcing the flavor off-diagonal axial 
transformations of flavor singlet operators, e.g. 

^r\P) = -^i^^^- These determine dp - bs, 
ds — bp, (It — bx, and rx- 

The remaining undetermined combinations are 
bs + bp + ds + dp, bT + dr and cIa (and ta)- 

We close by expanding upon one feature of our 
analysis. In the second step, when we enforce 
the PCAC relation, we need the 0{a) improved 
expression for quark masses 

(AM)=Z™ [{\M){\ + b^{M)) + b„,{\M^)] (11) 

[(M)(l + d„,{M)) + d^{M^)] . (12) 

Note that singlet and non-singlct masses are nor- 
malized differently, so that individual masses have 
off-diagonal renormalization even at 0(1): 

m, = Z„ K- + (r,„ - 1)(M)/A^/ + 0{a)] . (13) 

This effect is absent in the quenched approxi- 
mation, where = 1, as well as c?„i = bm, 

bm — dm — 0. 



*Some of combinations of coefficients depend on Nf, and 
we quote results here for Nj = 3. 



It turns out that these new improvement co- 
efficients are related to those we have previously 
introduced, and we need to determine these rela- 
tionships in order to extract all relevant informa- 
tion in step 2 above. In other words, we need to 
generalize the quenched relations Z^Zs = 1 and 
bs = —2bm @- To do so we note that it is consis- 
tent with the WI to take the matrix elements of 
the improved renormalized scalar density to be 



dniH 



dfhi 



(14) 



with H an arbitrary hadronic state. Note that 
the lattice spacing, and not is held fixed in 
the derivative. The r.h.s. can be evaluated in 
terms of derivatives known from the form of the 
lattice action 



dniH 



^druH 



dgo 



(H\S^"^\H), (15) 



{H\U{F^,F^,)\H). (16) 



This results in the relations* ZsZ„i = 1, rsr„i = 
1J>S = -2bm, bs ^-b,n, 3/s = 2{b„i - dm), and 
3ds = 36m + 26m — 6rf,„ — 2d,„. One also finds two 
constraints, namely ds = 36^ + 65, and 255 = 
bg/gQ. The latter can be used as a consistency 
check on the calculations of gs and bg. 
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